LAB THREE: SELECTIVE EXCITATION¶

This lab covers selective excitation, RF pulse design using Sinc and SLR, slice profiles, slice offsetting, and the relationships between pulse bandwidth, gradient strength, slice thickness and SNR.


Required Resources¶

For this lab you will need:

  • ilumr fitted with 10mm RF Probe
  • Shim Sample (#001)
  • Phantom #005 (Three Cell Phantom)
  • Phantom #006


Setup Task: Run the Notebook¶

  1. Edit the cell below to set the LAB_USER_NAME variable to your name
  2. Click Run->Run All Cells in the in top menu bar of jupyterlab
  3. Open the Table of Contents side-bar on the left edge of jupyterlab to aid in navigation

Important: To initialise this notebook, edit the cell above to set LAB_USER_NAME to your name, then click Run->Run All Cells in the top menu bar.

Data will be saved to /home/data/REPLACE_ME/Lab 3

1. Projection Images¶

Thus far, the 2D images created in the labs have been projection images. These images are created using a rectangular RF-pulse (known as a hard pulse) and both frequency and phase encoding gradients (Fig.1.1). The frequency and phase encoding gradients allow the aquired signals to be sorted into 2D space to produce an image. The orientation of the image will depend on which axes the frequency and phase gradients were applied along. Since there is no form of encoding in the third dimension of the sample, the signals along this axis will be summed in the 2D image.

No description has been provided for this image
Figure 1.1: Setting the frequency and phase encoding in the x and y axes produces a projection image where the varying signal components in the z direction are summed into average values.

Task 1.1: Generating a Projection Image¶

  1. Insert Phantom #005 at the correct depth.
  2. Choose XY for the phase encoding and frequency encoding axes and start the experiment by pressing "Run". (If the image is offset, run the Find Frequency dashboard)
  3. Question: Can you tell the vertical position of the cells?
  4. Try changing one of the encoding axes to Z to get a side profile. If the sample is off-centre, use this image to adjust it.
  5. Insert Phantom #006 and repeat the XY and ZY images.
  6. Question: Is the text embedded in the phantom visible/readable?

As you can see, projection images aren't always useful because summing a large amount of signal in the object's third dimension makes it almost impossible to see details that are only present in specific parts of the sample. If we want to see detail in a specific portion of the object we need to be able to limit the amount of signal that is being summed in one direction. This is done through slice selection.

2. Frequency Selective Excitation¶

In order to image a single slice rather than the entire sample, two main changes need to be made to the pulse sequence:

  1. A slice-select gradient must be added along the axis perpendicular to the plane of the desired slice. This causes the spins along the chosen axis to precess at different frequencies making it possible to selectively excite them by transmitting at a specific frequency.

  2. The RF excitation pulse must be modified to only transmit a small range of frequencies. So far, we have used short (~20 us) high amplitude pulses which excite a very broad range of frequencies (called a "hard pulse"). As the pulse width is increased, the range of frequencies contained in it (the bandwidth) decreases, so we need to increase the pulse length to achieve a narrow slice. Flip angle is proportional to the area of the pulse, so a long, narrow-bandwidth pulse needs less amplitude to achieve the same flip angle. This is called a "soft pulse".

The resulting pulse sequence is shown in Figure 2.1 below.

No description has been provided for this image
Figure 2.1: Slicing along the z axis allows the sample to be imaged at a specific cross section index

The relationship of the slice thickness ($\Delta$z) to bandwidth ($\Delta$F) and slice-select gradient strength ($G_{ss}$) is governed by the following equation:

$$\Delta z = \frac{\Delta F}{\bar\gamma \cdot G_{ss}} \tag{1}$$

where $\bar\gamma$ is the gyromagnetic ratio (for hydrogen nucleus $\bar\gamma = 42.58 \times 10^{6} \ \mathrm{Hz}/\mathrm{T}$). Figure 2.2 shows examples of how slice thickness can be controlled by varying $\Delta F$ or $G_{ss}$.

No description has been provided for this image
Figure 2.2: How the slice thickness is related to the RF pulse bandwidth and slice-select gradient strength

In the next experiment, we will use a simplified pulse sequence (Fig.2.3 below) that performs a spin echo with a constant gradient present during both the RF excitation pulse and the readout. This allows the sample to be selectively excited and then imaged along the same axis to directly measure the excitation profile.

No description has been provided for this image
Figure 2.3: Pulse Sequence for the Constant Gradient Spin Echo experiment

Task 2.1:¶

  1. Insert the Shim Sample (#001) and run the below experiment with a pulse width of 20 us. You will get a 1D Z axis projection image of the sample.
  2. Increase the pulse width and rerun the experiment. Try with various pulse widths up to 1000 us.
  3. Note: With a short RF pulse, there should be roughly equal signal from all parts of the sample along the Z axis.
  4. Question: With a long RF pulse, do you still get equal signal from the whole sample? If not, which region dominates?
  5. Question: How does the Z thickness of the dominant region (slice thickness) vary with the pulse width?


Task 2.2: 2D RARE Using Rectangular Soft Pulse¶

  1. Insert Phantom #005, and run the experiment with the pulse width set to 20 us and the axes set to X-Y to get a projection image.
  2. Try increasing the pulse width to selectively excite only the middle cell. You can also switch one of the axes to Z to visualise the excitation profile as in the previous experiment.
  3. Question: Approximately how long does the RF pulse need to be to only get noticable signal from one cell?
  4. Question: What is the thickness of this slice? (You can check with the Constant Gradient Spin Echo experiment above)
  5. Question: If you continue increasing the pulse width, how is the SNR of the image related to the slice thickness?

3. Soft Pulse Design¶

In practice, a slice will have some finite thickness. Choosing a slice thickness is a trade off between two main factors: A thicker slice will have more signal because there are more nuclear spins available for excitation; a thinner slice will have less averaging over the slice axis resulting in better selection of individual features in the sample. In the previous experiments, a rectangular RF pulse with a constant amplitude was used. This resulted in a Sinc function shaped excitation profile (or slice profile). The side-lobes of the slice profile partially excite the top and bottom cells causing them to show faintly in the image.

For small tip angles ($\alpha \ll 90^\circ$), the small-tip-angle approximation states that the signal spectrum will be approximately proportional to the spectrum of the RF pulse, and therefore will have the same shape. If the tip angle at any point in the sample approaches or exceeds $90^\circ$ this approximation breaks down and the signal spectrum will be distorted, compared to the RF pulse spectrum, due to the sine relationship between signal and tip angle (Fig. 3.1).

No description has been provided for this image
Figure 3.1: How the RF pulse shape relates to the NMR signal

To perfectly isolate the middle cell, it should be fully excited (90 degree tip angle throughout) while the outer cells are not excited at all (0 degree tip angle). This requires the frequency spectrum of the RF pulse to have a rectangular shape, which can be theoretically achieved with a Sinc shaped RF pulse amplitude modulation as shown in Figure 3.2:

No description has been provided for this image
Figure 3.2: An ideal RF pulse creates a rectangular slice select profile

In practice, the theoretical Sinc function is infinitely long, and must be truncated to be used. This may be done at the zero crossings to minimise ringing artifacts (Gibbs phenomenon). Additionally, the Sinc shape can be multiplied with a gaussian to further reduce the Gibbs phenomenon. This process is called Apodization.


Task 3.1: Visualising Sinc Pulse Design¶

The experiment below is designed to help you visualise how changing the parameters of a soft pulse changes the profile in the frequency domain. The RF pulse you design will be used for the following experiments.

  1. Experiment with the three soft pulse parameters and observe how the shape of the RF-pulse and slice profile changes.
  2. Question: How does changing the lobe width affect the spectral bandwidth?
  3. Question: How does adding more lobes improve the shape of the slice profile? What is the trade-off regarding the overall pulse width?
  4. Question: What happens when apodization is too large? Pick a value that looks like a good compromise.


Task 3.2: 1D Slice Profile¶

With the Sinc pulse designed, we can now use it to create a 1D image of the slice profile. You will need Equation 1: $$\Delta z = \frac{\Delta F}{\bar\gamma \cdot G_{ss}}$$ where $\bar\gamma = 42.58 \times 10^{6} \ \mathrm{Hz}/\mathrm{T}$ for the hydrogen nucleus.

  1. Question: Based on the bandwidth ($\Delta$F) of the Sinc pulse you have designed and a slice select gradient $G_{ss}=100 \ \mathrm{mT}/\mathrm{m}$, calculate the expected thickness of the slice profile.
  2. Insert the Shim Sample (#001) (make sure to centre the sample using the depth gauge), set the slice select gradient, and run the experiment.
  3. Question: How does the width of the slice profile compare to your predicted width?
  4. Try changing the flip angle to 30 and 180 and compare the resulting spectrum.
  5. Question: At what flip angle (out of 30, 90, or 180) does the signal spectrum look most similar to the RF pulse spectrum shown above?

Note: Phase Rewinding Gradients¶

There may be some twisting in the imaginary components of the spectrum in the previous experiment. This is due to dephasing. The application of the slice-select gradient during the soft pulse causes the spins to accumulate different amounts of phase based on their position along the the gradient's axis. This phase dispersion of the transverse magnetization will result in a loss of signal when averaged together. Since the dephasing that occurs is linear, the phase effects can be cancelled by applying a gradient in the opposite direction. This is done using a rephasing gradient. In general, the area of the rephasing gradient is half that of the slice-select gradient. This is done under the assumption that the majority of the spins are tipped into the transverse plane at the centre of the 90-degree pulse and, therefore, phase dispersion primarily occurs during the second half of the slice select gradient. While this assumption is a good starting point, the size of the refocusing gradient needs to be calculated based on the soft pulse design if the phase effects are to be cancelled perfectly.

No description has been provided for this image
Figure 3.3: Slice Select Spin Echo pulse sequence without (left) and with (right) rephasing lobes

Task 3.3: 2D RARE Using Sinc Shaped Soft Pulse¶

  1. Insert Phantom #005, and run the experiment below with the Z-Y axes selected.
  2. Adjust the slice select gradient to change the slice thickness until only the middle cell is visible. You may need to also modify the sinc pulse design above if the gradient is reaching the maximum allowed. The sample can also be moved up or down if it is not centred.
  3. Change the axes to X-Y and run again. You should see only the middle cell appearing in the image. If not, go back and adjust the soft pulse/gradient parameters until you do.
  4. Save/screenshot the image

4. SLR Soft Pulse Design¶

The Sinc method of designing a slice select RF pulse works well at small flip angles, but is not ideal for 90 degree pulses and is quite poor for 180 pulses. A more optimal method, that is widely used, is the Shinnar-Le Roux algorithm. This algorithm takes into account the non-linearity of NMR excitation to achieve equal ripple in the final magnetization spectrum and, therefore, signal spectrum.

For 90 degree pulses, this method results in a RF pulse spectrum with minimal sidelobes at the expense of poor flatness on top. However, since the signal level resulting from a near 90 degree flip angle pulse is not very sensitive to the exact flip angle, the slice profile is still quite flat.

The SLR algorithm is also used to design asymmetrical RF pulses which are used for advanced techniques.


Task 4.1: Pulse Design For a Target Slice Thickness Using SLR¶

You will need Equation 1: $$\Delta z = \frac{\Delta F}{\bar\gamma \cdot G_{ss}}$$ where $\bar\gamma = 42.58 \times 10^{6} \ \mathrm{Hz}/\mathrm{T}$ for the hydrogen nucleus.

  1. Question: Calculate the required bandwidth to achieve a 1 mm slice with the maximum slice gradient, $G_{ss}=150 \ \mathrm{mT}/\mathrm{m}$, and enter the value below (take care with units).
  2. Insert the Shim Sample.
  3. Measure the slice profile with the above slice select gradient and check the thickness is correct. You may increase the pulse duration to improve the sharpness of the edges of the slice.

5. Controlling Slice Position¶


Task 5.1: Locating the Hidden Text¶

  1. Insert Phantom #006 and adjust the settings to get an XY image of a slice through the text inside the phantom.
  2. Hint: Try reducing the slice select gradient to get a thicker slice and viewing a Z-Y projection to identify the plane where the text is. You can measure the Z offset of a bright plane and calculate the required frequency offset to position the slice correctly. The text is located in one of the bright planes that are visible in a Z-Y projection.
  3. Save/screenshot the image